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Abstract 

We show (in a completely analytical and exact manner) that an efficient magnetic field ampli- 
fication method is operative during the bounce in a time-dependent gauge coupling model. The 
cosmological magnetic fields so generated have particular spectral features, and may be observed 
by future CMB measurements and by direct cluster measurements. 
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Introduction 



The origin, evolution, and structure of large-scale magnetic fields are among the most 
important open issues in astrophysics and cosmology. The magnetic field present in galaxies 
is typically of the order of a few /iG, and is coherent on the galactic scale. In the case of 
clusters, the observed field is of the order of the fiG and is correlated over 10-100 kiloparsecs. 
The standard mechanism to account for these fields is the dynamo which amplifies small 
seed fields to the abovementioned values [2j. Three different origins for the seeds that 
fuel the dynamo have been discussed. They can be called astrophysical, cosmological, and 
primordial. The most popular astrophysical mechanism for the generation of these pre- 
galactic magnetic fields is the Biermann battery [sj, based on the conversion of the kinetic 
energy of the turbulent motion of the conductive interstellar medium into magnetic energy. 
The battery works for the generation of seeds to be amplified by the dynamo in galaxies, 
but can hardly account for the fields in clusters [4, 0]. Other astrophysical mechanisms 
(involving for instance starbursts, jet-lobe radio sources, or accretion disks of black holes in 
AGNs seem to need pre-existing magnetic fields. 

The cosmological mechanism is based on the generation of magnetic fields from cosmo- 
logical perturbations The idea is that due to collision effects, there is a difference in 
the velocity of protons and electrons around and after the decoupling time, which leads to 
an electric current that generates a magnetic field. By this mechanism, magnetic fields of 
« 10~ 18 G at 1 Mpc scale and « 10~ 14 G at 10 kPc scale can be generated fl. 

The current prevalent view is that the magnetic fields observed in galaxies and clusters 

n 

have a primordial origin (see [10J for a review). The processes that may account for this 
origin can be divided in two types: causal (those in which the seeds are produced at a given 
time inside the horizon, like QCD and EW phase transitions [k]]), and inflationary (where 
correlations are produced outside the horizon ll|). The former presents problems due to the 



fact that the coherence length of the fields so generated is much smaller than what is needed 
as input seed fields for a galactic dynamo [6] , while in the latter, vacuum fluctuations of the 
electromagnetic field are "stretched" by the evolution of the background geometry to super- 
horizon scales, and they could appear today as large-scale magnetic fields. However, since 
Maxwell's equations are conformally invariant in the FRW background, the amplification 
of the vacuum fluctuations (which amounts to particle production) via inflation can work 
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only if conformal invariance is broken at some stage of the evolution of the universe. There 
are several ways in which this invariance can be broken Il2|: non-minimal coupling between 

nnn 

gravitation and the electromagnetic field 1131 . 114J ] , quantum anomaly of the trace of the 
stress-en ergy tensor of electrodynamics [15J, coupling of the EM field to a charged scalar 



field fill . Il6| |. exponential coupling between a scalar (whose potential drives inflation) and 



the EM field 



171 ] , and a non-zero mass for the photon [18( . Yet another possibility is dilaton 
electrodynamics [19], in which there is a scalar field (the dilaton) which couples exponentially 
to an Abelian gauge field. In this model the inflationary expansion is driven not by the 
potential but by the kinetic term of the scalar field. The exponential coupling is naturally 
implemented in the low-energy limit in string theory [20] , and in Weyl integrable spacetime 
(WIST) [21], and can be viewed also as a time dependence of the coupling constant, an 
idea considered first by Dirac [22]. This avenue has been pursued by Giovannini in a series 



of articles 



231 ] . More recently, a model with both a dilaton and an inflaton was analyzed 



25J. In a 



1 these articles, 



261 ] the features 



in [24j], and later generalized to a noncommutative spacetime 

different aspects of photon production have been analyzed. We studied in 

of the creation of photons in a non-singular exact solution of Einstein's gravity coupled to 

the dilaton, in which the passage through the bounce is described by the equations of the 

model without resorting to unknown (Planck- scale) physics [36]. The conformal symmetry 

was broken through the exponential coupling of the dilaton to the EM field. Using the 

I I 

formalism of the squeezed states [27], the number of created photons in this model was 
obtained through numerical calculation in 26]. In this article we we shall re-address the 
creation of photons in the above mentioned model by a different quantization technique 
(i.e. canonical quantization). The advantage in this case is that the results are obtained in a 
completely analytic manner. Moreover, we present here a complete analysis of the produced 
magnetic field spectrum. We shall find an exact solution for the mode equation of the vector 
potential A^, and from it an exact expression for the spectrum of the magnetic field will 
follow. We shall see that the requirement that the magnitude of the seed fields produced by 
the bounce are compatible with observation can be fulfilled by a large range of values of the 
parameters of the model. 
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I. FIELD EQUATIONS 



A time-dependent gauge coupling is a generic feature of 4-dimensional theories obtained 



by compactification of a more general theory such as Kaluza-Klein theory_ 
theory 



281 ] and string 



29| | . Time- dependent gauge couplings are also present in WIST [21J. In all these 
cases, the action can be conveniently written in the form 

S=\f d'x^Tg f(oo)F aP F aP , (1) 

where uo is the the dilaton (in the case of string theory) or the scalar field associated to Weyl 
geometry (in the case of WIST) and F a p is an Abelian field. The function f(u>) will be set 
in the following to e~ 2uJ , which corresponds to the case of string theory and WIST. 
From now on we shall work in conformal time, with the FRW metric given by 

ds 2 = a(r]) 2 [dt] 2 — r )ij(x)dx i dx : ' 1 \ , 

where 7^ is the metric of the hypersurface 77 = constant, and i — 1, 2, 3. 

Using the definition of the canonically normalized vector potential A^ = e^A^, and the 
radiation gauge, defined by A$ = 0, V.A = 0, the Lagrangian density in terms of A a can 
be obtained from Eqn.([T]): 

C = \[{A!tf + {u' 2 -u")AZ-{d t A a ) 2 ], 

where the prime indicates derivative w.r.t. conformal time rj. The equation of motion that 
follows from this Lagrangian is 

uAi + (w" - uj ,2 )Ai = 0. (2) 

Since we are dealing with vector quantities we shall use a vector basis P l defined on a 
hypersurface i] = constant by the following relations [3^] 

p i = P^x*), Y j V< Vj P l = -m 2 P l , Y j ViPj = 0. 

The eigenvalue m denotes the wave number of the corresponding vector eigenfunction of the 
Laplacian operator on the spatially homogeneous hypersurfaces. The spectrum of eigenval- 
ues depends on the 3-curvature e, and is given by 

m 2 = s 2 + 2, < s < 00, for e = -1, 



m = s, < s < oo, for e = 0, 
m 2 = s 2 - 2, s = 2,3,.., for e = 1. 
In terms of the basis P % the vector potential can be expanded as follows: 

^(r / ,x) = (2vr)- 3 ^ J d 3 m A^M 

This expression is valid in the cases e = —1,0. For e = 1, the integral must be replaced 
by an sum. The index a = 1,2 describes the two transverse degrees of freedom. We shall 
work with an expansion in terms of travelling waves, hence a takes the values "+" or "— " 
according to 



ma \ i 



In the case of a 3-space with zero curvature, the P % are such that P^ l {x) = e l a e ±m ' x , with 
e a .m = 0, e a .e a> = 8 aa i. Substituting this expansion in Eqn.([2]) we obtain the equation 
that governs the evolution of each of the polarization modes of A 1 : 

A^\v) + [m 2 -V( V )]A^( V ) = 0, (3) 

and a similar equation for Aml(rj), where V{rj) = —uj" + u' 2 . 

The vector potential can be quantized following standard procedures, by transforming 
the mode functions into operators [31] : 



Here 

n m (v) = •KSn) , [«« , akJ = S aP S {3) (m - p) , 
and the modes A m obey Eqn.(j3]). 

II. THE BACKGROUND 



As discussed in [32| , there are nonsingular solutions in the theory of WIST that describe 
a FRW geometry plus a scalar field, as well as in string theory [33]. Let us briefly review 



these solutions. The EOM for gravitation plus scalar field written in conformal time are 32] 



A 2 

a' 2 + ea 2 H (J a) 2 

6 v ' 



0. 



to 



7« 2 , 



(4) 



where 7 = constant, and A 2 is the coupling constant of the scalar field to gravity. Notice 
that the equation of state of the scalar field is given by p w = p^, where 



Pu 



A 2 fu' 



From Eqns.(jl]) we get 



,'2 



-ea 



(5) 



where we have defined a 2 , = X'y/y/Q. Eqn.Q shows that only solutions with e = — 1 are 
possible. Hence, from now on we shall restrict to the negative curvature case, for which 

n 

Eqn.([5]) can be easily integrated. The result of the integration for the scale factor is |37| 



Eqn.fll) yields 
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a{j]) = ao^/cosh.(2r] + 5). 



(7 ? ) = ±^arctan(e 2 ^ 5 )+4^- 
KU 2A V 1 v^A' 



(6) 



(7) 



where 5 is an integration constant. Notice that the model here presented has two independent 
constants. The plots for these functions are given in FigJU The scale factor displays a 
bounce, produced by the violation of the strong energy condition by the scalar field 




FIG. 1: Plot of the scale factor and u as a function of the conformal time for A = 00 = 1. 
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III. MODE EQUATION AND SPECTRUM 



In terms of the field given in Eqn.Q, the potential takes the form 

2crsinh(2r/) + a 2 



v(v) 



cosh 2 (277) 



(where a = V6/X) and its plot is given in Fig. ([2]). In order to calculate quantities of interest 




FIG. 2: Plot of V(rj) for A = 0.1, 0.3, 0.5. 



related to the magnetic field, we need to solve the mode Eqn.Q, which in the case at hand 



is 



AJjlf +[m 2 - 



2Av / 6sinh(2r/) W 



A 2 cosh 2 (2r?) 

The general solution of this equation is given by 



A m (v) = 0. 



Am(r]) = diZ^rj) + d 2 Z 2 (r]), 



where 



Z 1 ( V ) = e-^F(a, -a;c;y( V )), 
Z 2 (v) = e^yivY^Hl + a-c,l-a-c;2-c; 2/(77)), 
with oj{rj) given by Eqn.Q and 



ira 



1 + ia 



y 



l + i sinh(2r/) 



2 7 2 1 2 

The coefficients in the superposition in A m {v) have been determined in order to satisfy the 
following asymptotic conditions: 

AJji) = -== e~ imr] = Ai n \r]) for t] -> -00 



'2m 



7 



(which implies that |A m (?7)| 2 = and 



They are given by 



[r(i= 



-imr) 



1—im—ia \ p/ l — im+itr N 



'2m 



sinh( 



7T(7 
2 ' 



cosh(^f) 



imri 



'2m 



= At ut \v) for r]^oo. 



H^^f. T ^ + irn+ia ^ p^ l+im-ia ) r ( l-im+i<0 



do 



m r(i^) r(^p) [ r( 1+ y g ) r( 1 -^ g ) + r( 1 ~^ +i(T ) r( 1+i ™~ ia )] ' 

-vWe^T+f) r( 1+ y ig ) r(i±np^) v( l - ir ^ ia ) 
V2 r(^) r(i^) [ r(i±i22±i£) r(±^p^) + r(i±np^) f(^f±^j] ' 



From them we can read out the Bogoliubov coefficients that relate the ingoing with the 
outgoing fields, respectively given by 

A^\x, v ) = J d s mX)[^.^ ) (»7)^. + a^^*(i7)*P ia ,J, 



J s 

The operators a^ s and are related by the Bogolubov coefficients and (3^ as follows: 

In our case matrices aij and /3y are given by 

*jrn = S(j-m)ct, (3^ = 5(j+m)(-iy C -\ 

where 



[r(i^)j- 



-"m p / 1— im— to- \ p / 1— im+icr \ ' 



"(-I) 



2 ,n 2 

sinh(7T(j/2) 



(8) 
(9) 



with lei; 1 2 - 1 2 



cosh(7rm/2) 

r-,,,1 - -= 1. These coefficients yield the spectrum of the magnetic field through 



the expression 



hm 4 c 

pym) 



a 4 7r 2 







lc"| 2 - 




1 L 'ml 





„+ -2imr) , *+ 2imr) 



which can be rewritten as 



JytTl C 

= TI-T (1 + 2| c~ | 2 - 2|c~| |c+| (cos(2m?? - m ))), 



a 4 vr 2 



8 



where 9 m = arctaii (S(c+)/5i(c+)). When na 3> 1 (that is, when 3> A), and using 

l c m| 2 ~ \ c m\ 2 = 1- ^is expression becomes: 



plus corrections 0(l/\ c _ | 2 ). This expression for the spectrum is valid for i] ^> 1. When 
translated into cosmological time, it gives t ^> a x a , where a ~ 10 2 . Consequently, our 
result is valid for t = today if we choose a convenient value of a . 

It follows from the spectrum that the amplification factor with respect to the conformal 
vacuum peaks for the modes with momenta such that m ~ 1.31, and it is of the order, 



which is exponentially large for A< 1. Using the constants ao and A (linked to 7 through 



where t r is the time at which the scalar field is negligible, in such a way that the electro- 
magnetic field is free again. 

At a comoving scale of about 10 kpc the strength of conformal vacuum fluctuations is 
of the order 10~ 55 gauss. To get it to the strength required to seed the galactic dynamo, 
.Bseed ~ 10~ 20 gauss, which is a conservative estimate, we get for the required amplification 
factor, exp(7r^6/A) ~ 10 35 , or A ~ 0.1. If one takes for the comovine; scale the size of the 
Universe, 4 x 10 6 kpc, the amplification factor becomes, exp(7TA/6/A) ~ 10 46 for which we 
find A ~ 0.07. 

IV. DISCUSSION 

We have shown that when the universe evolves through a bounce in a theory in which the 
EM field is coupled to a scalar field through an exponential (as is the case of string theory 
and WIST), EM fields of the magnitude needed for a seed can be generated for a large range 
of values of the parameters in the theory (A and ao). The expression for the spectrum has 
been obtained un a complete analytical manner. Notice that the evolution of the geometry 




(10) 




(11) 
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on this model is continuous, since the bounce is caused not by "new physics" but by the 
scalar field. 

As regards the matching of the bounce era onto radiation and matter era, since gauge 
fields are conformal and so is the bounce at asymptotically late times, the matching onto 
the radiation era solutions, 

' :e +mv and -^Le^"" 7 (12) 



V2m V2m 

is trivial (there is no Bogoliubov mixing) , which in turn implies that the spectrum calculated 
during the bounce remains unchanged during radiation and matter eras. 

The model presented here shows that the mechanism for production of magnetic fields 
during the bounce may work, but we do no claim that it actually solves the problem. Some 
more work is needed in certain aspects. In particular, if we want to interpret the spectrum 
obtained here as that observed today, we need to assume that no further amplification occurs 
at the bounce-radiation transition. The incorporation of matter into the model should also 



be considered (perhaps along the lines of [26|. We hope to deal with these matters in a 
future publication. 
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